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We present an investigation of the total decay rate of the (ground state) B c meson within 
the framework of the relativistic constituent quark model formulated on the light-front 
(LF). The exclusive semileptonic (SL) and nonleptonic (NL) beauty and charm decays of 
the B c meson are described through vector and axial hadronic form factors, which are 
calculated in terms of quark model LF wave functions. The latter ones are derived via 
the Hamiltonian LF formalism using as input the update constituent quark models. The 
inclusive SL and NL partial rates are calculated within a convolution approach inspired by 
the partonic model and involving the same B c wave function which is used for evaluation 
of the exclusive modes. The framework incorporates systematically 84 exclusive and 44 
inclusive partial rates corresponding to the underlying b — > c and c — > s quark decays. We 
find tb c = 0.59±0.06 ps where the theoretical uncertainty is dominated by the uncertainty 
in the choice of LF wave functions and the threshold values for the hadron continuum. 
For the branching fractions of the B~£ — * J/i/jfj^Un and B£ — * J/tpir + decays we obtain 
1.6% and 0.1%, respectively. 
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The theoretical interest in the study of the B c meson, the bound state of the be system with 
open charm and beauty, is stimulated by the experimental search at CDF and LHC Recently 
the CDF Collaboration reported the observation of B c in 1.8 TeV pp collisions at Fermilab JT|. 
The CDF results for the B c mass and lifetime are M Bc = 6.40 ± 0.39(stat) ± 0.13(syst) GeV/c 2 
and t Bc = 0A6t° il(stat) ± 0.003(syst) ps. The physics of B c mesons has stimulated much 
recent works on their properties, weak decays and production cross section on high energy 
colliders. For a review see . 

Similar to D and B mesons the ground be state is stable against strong or electromagnetic 
decay and disintegrates only via weak interactions. The weak B c decays occur mainly through 
the CKM favored b — > cW + transitions with c being a spectator, leading to final states like 
J/i/jiu, and c — > sW + transitions with b being a spectator, leading to final states like B s 7r, 
B s iv. Weak decays of charmed and bottom hadrons are particularly simple in the limit of 
infinite heavy quark mass, where the decay rate of a hadron Hq containing a heavy quark Q 
is completely determined by the decay rate of the heavy quark itself. In this limit, one might 
expect that T(B C ) « T(B°) + T(D°) yielding tb c ~ 0.3 ps, with c-decay dominating over b- 
decay In reality, hadrons are bound states of heavy quarks with light constituents. The 
account of the soft degrees of freedom generates important pre-asymptotic contributions due 
to the Fermi motion of a heavy quark inside the hadron. These effects have a significant impact 
on the lifetime and various branching fractions of B c . The various calculations of tb c have been 
reported in the literature @-||. The wide range of predicted lifetimes r Bc = 0.4 — 0.9 ps, 
reflects the uncertainty due to the various model assumptions on the modification of the free 
decay rates due to the bound state effects and the limited knowledge of the heavy quark masses. 

In this paper, we use an approach in which non-perturbative QCD effects are mocked up 
by a light-front (LF) wave function of the hadron ||. The internal motion of a heavy Q— 
quark inside the heavy flavour meson Hq is described by the distribution function F(x) = 
J d 2 p±\t/j(x, p\)\ 2 , where |-^(x,p5_)| 2 represents the probability to find a quark Q carrying a LF 
fraction x = Pq/ P Bq of the meson momentum and a transverse relative momentum squared p\ . 
A relevant feature of this approach is that both exclusive and inclusive decays are coherently 
treated in terms of the same heavy quark wave function. So far the approach has been applied 
only for the exclusive and inclusive partial widths of B° |3[], where it has been found that the 
overall picture is quantitatively satisfactory. Here we extend previous calculations started in 
Ref. (?j to compute the lifetime and various decay branching fractions of B*. Our aim is to 
constrain the model dependence of the calculated tb c related to various choices of ^(x,pjj. To 
this end we will make use of different LF wave functions, constructed via the Hamiltonian LF 
formalism (see e.g. ||) adopting recently developed relativized and a non-relativistic quark 
models |,[1|. 

To start with, we briefly remind the main points of the procedure used to calculate heavy 
meson partial widths. For more details see Refs. |J, 0. Consider the SL decay rates first. 
Instead of considering the exclusive modes individually we will sum over all possible hadronic 
final states X. This sum includes hadronic states with a large range of invariant mass Mx- 
For the heavy-flavour mesons, the energy which flows into hadronic system is typically much 
larger than the energy scale Aqqd which characterizes the strong interactions. Consequently, 
for a wide range of the phase space an inclusive description based on quark-hadron duality is 
appropriate. It will be valid over almost all of the Dalitz plot, failing only in the narrow corner 
region where the observed mass spectrum is dominated by the two narrow peaks corresponding 
to the transitions Hq — > P and Hq — > V, where P and V are the lowest lying pseudoscalar 
and vector mesons. Accordingly, the total SL rate of the Hq meson has been represented in 
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the following hybrid form 



r(H Q -> xtvt) = t(h q -> p^) + r(# Q -> + v(H Q -> x'^), (1) 

where X' represents the hadron continuum including also the resonance states higher than P 
and V. The usefulness of such an expansion rests on large energy release in the inclusive decay. 

Our calculations of the exclusive rates T(Hq — > Plvi) and T(Hq — > V^) use the hadronic 
form factors that depends explicitly on dynamics of specific channels. The relevant formulae 
valid also in the case when the lepton masses are not negligible are collected in |7| and we 
do not quote them here. Instead, we will concentrate on the calculation of the inclusive rate 
T(Hq — > X'tvf). The modulus squared of the amplitude summed over the final hadronic states 
is written as \M\ 2 = (G 2 F /2)\VQiQ\ 2 L a(3 W al 3, where Vq<q is the relevant CKM matrix element, 
L a P is the leptonic tensor and W a p is the hadronic tensor. We use the notation of Ref . [ I 1 



The hadronic tensor can be expressed in terms of five structure functions W\ to W5 that depend 
on two invariants, q 2 and q$, where q is the 4-momentum of a lepton pair and go is related to 
Mx by: go = {Mr + <? 2 ~ M x )/2Mh q . It is convenient to scale all momenta by Mh q , letting 
M\ = s ■ Mjj Q and q 2 = t ■ M^. Evaluating the contraction L a ^W a p we arrive at the formula 
for the SL inclusive width 

TsL = ^r \V Q/Q \ 2 J S L I dt®(t,ml,ml) f ds-^-G(t,s), (2) 
6 J J M Hq 

where the prefactor T = (G'pMfj )/(47r) 3 sets the overall scale of the rate, and Jsl ~ 0.9 rep- 



resents the effect of the radiative corrections |I2" |. In Eq. (0) $(g 2 , mf, mf) = J 1 — 2A + + A 2 , 
with A± = (m 2 ± ml)/q 2 , being the lepton masses, Ai = A + — 2A 2 , A2 = A + — A 2 . Fur- 



thermore, 2|q|/M# Q = a(t, s) = y(l + t — s) 2 — At, and the limits of integrations in the t — s 

plane are given by s min = (M^/M Hq ) 2 , s max = 1 - </t, t min = (m 1 + m 2 ) 2 /m 2 Q} t max = 

(1 — M x /Mh q ) 2 , where M x is the threshold value at which the hadronic continuum starts. 
The function G(s, t) is expressed in terms of the linear combination of Wi(t, : 

G(t, s) = 3t(l + Ai - 2\ 2 )W 1 + ( (1 + Ai)^- + |a 2 ^ W 2 + ^\ 2 t ((I + 1 - s)W 4 + tW 5 ) . (3) 

To calculate the structure functions Wj we use the parton approach of Refs. [|| (see also |T3| ) 
based on the hypothesis of quark-hadron duality. This hypothesis assumes that the inclusive 
decay probability for which no reference to a particular hadronic state is needed equals to one 
into the free quarks. The basic ingredient is the expression for the hadronic tensor W a p which 
is given through the optical theorem by the imaginary part of the quark box diagram describing 
the forward scattering amplitude: 

f dx 

W af3 = J —w a/ 3(pQ,,pQ)6(eQi)5{(p Q - q) 2 - m 2 QI ]F(x), (4) 


where w a p(PQ',PQ) — § 2 UQ'O a Ub ■ uqO'Iuqi is the parton matrix element squared. In what 

spins 

follows we choose the purely logitudinal kinematics with qj_ = 0. In Eq. @ the 5-function 
1 The function W3 does not appear in the total SL width integrated over the lepton energy 
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corresponds to the decay of a Q-quark with momentum pq = xPh q to a Q'-quark with the 
momentum pg> = Pq — q and has two roots in x, viz. 5[(j>q — q) 2 — m 2 Q,} = [S(x — x + ) + 5(x — 
x-)]/(M% Q \x+-X-\), where x±(t,s) = |(1 + 1 - s ± ^(1 + t - s) 2 - At + Am 2 Q ,/M 2 HQ ) and the 

quark transverse momenta being neglected Q. The root x_ is related to the contribution of 
the Z-graph arising from the negative energy components of the Q'-quark propagator and is 
prohibited by the 6{eQi) in Eq. (U). Using the explicit expression for w a p in Eq. (|J) one obtains 

Wx = F(x+), W 2 = A * +F ^ X+ \ ) W 3 = W 4 = -2 l . , W 5 = 0. (5) 

— |— — %L/ _|_ iX- — 

In this way, the SL inclusive width gets related to the distribution F(x). Substituting Wi from 
(|5|) into (^) we obtain 

_ tmax Smax 

T(H Q -*X'£v e ) = -T J SL \V QQ/ \ 2 J dt$(t) J dsa(t,s) 

tmin s min 

[(1 + Ax)« 2 (t, a) - X + x + 3t(l - A 2 )}F{x+). (6) 

The calculation of the NL decay rate closely follows the SL one. We expect Hq decays 
to multimeson states to proceed predominantly via the formation of a quark-antiquark state, 
followed by the creation of the additional qq pairs from the vacuum. The effective weak La- 
grangian, e.g. for b — > cuq processes with q = d, s is given by = ^V c bV uq {ciOi + C2O2), 
where Oi and O2 denote current-current operators with the color-non-singlet and color-singlet 
structure, respectively. The lepton pair is substituted by a quark pair, and the Wilson coef- 
ficients Ci(/x) are the perturbative QCD corrections that describe the physics between the W 
boson mass and the characteristic hadronic scale \x ~ trq of the process. We shall use the 
values HH ci(mfc) = 1.132, 02(^5) = —0.286; ci(m c ) = 1.351, C2(m c ) = —0.631, obtained 



at next-to-leading order with the evolution of the running coupling constant being done at 
two-loop order using the normalization a s {m z ) = 0.118 ± 0.003. 

When calculating the NL decays we use again the hybrid approach. We write J dt in (0) 



imax 

as / dt = J dt5(t — £j>) + / dtS(t — Qy) + / dt, where t' min is now related to the threshold for 

f . 

mm 

the hadron continuum produced by the W current, and (py = M P y/M Hg , with M P y being 
the mass of a pseudoscalar or vector meson. The latter integral is treated in the same way 
as for SL decays. According to Eq. (P it is the sum of exclusive, rates and inclusive rates. 
The inclusive NL rate is given by Eq. @ with the substitution T Jsl — > 3r |V^ g | 2 ?7 where 
rj = C 2 + (N C + 1)/(2JV C ) + c 2 _(N c - 1)/(2N C ), and c± = a± c 2 . In the limit of large number of 
colors N c which we use below 77 = c 2 + c\. 

Additional contributions arise from the first two terms containing S(s — Cpv)- ^ OT ^ ne 
exclusive two-meson decays Hq — > PP, PV, VP, and VV we use the factorization approach 



15| with the effective QCD coefficients a\ = 1, a 2 = —0.3 and a\ = ci(m c ), a 2 = C2(m c ) for 



the b c and c — > s transitions, respectively 0. The partial width of the inclusive NL decay 

2 By the quark masses toq and mg/ we hereafter understand the "constituent" quark masses taken from a 
particular constituent quark model. 

3 The coefficients 011,2 correspond to two flavour-flow topologies relevant for our discussion: the so-called 
"tree topology" (class I amplitudes), dominated by color allowed external W-dec&y and the "colour-suppressed 
tree topology" (class II amplitudes), dominated by colour-supressed internal W-decays There is the third 
class decays in which the a± and amplitudes interfere. These decays play the important role in the case of 
D mesons, but are less important for B c and are absent for B°. 
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Hq — > P + Xq,q in which the final state contains a charged or neutral pseudoscalar meson 
directly generated by a color-singlet current are given by 



(i-C) 2 

r ^ = ^72 — T o\V qiq2 \ \V Qi q 2 \ ■ / + - (x_ +3a; + )C ]— p ^-j— . (7) 

Hq , J \ X + — X -\ 

The analogous inclusive width Tyx for the production of a vector meson is obtained from ([/]) by 
the substitution fp — > fy and [x + (x + + x_) 2 — (aL + 3x + )( 2 ] — ► + a;_) 2 — (3aL + a: + )£ 2 ]. 

Here /^y are the pseudoscalar and vector meson coupling constants to the W current, and 
x± = x±(( 2 , s). The constants fpy are taken from Ref. ||16| . For the B c meson mass we use 
Mp c = 6.3 GeV, the other meson masses are taken from the recent PDG publication |17 |. 



The non~perturbative ingredient of our calculations is the LF wave function ijj(x,p\). In 
what follows, we will adopt for the latter the functions corresponding to the various equal time 
(ET) quark model wave functions. There is a simple operational connection between ET and 
LF wave functions ||, which allows to convert the ET wave function w(p 2 ) into relativistic 
LF wave function ^(a^pjj = ^p, where p 2 = p]_ + p 2 -, p z is the longitudinal momentum 



/4i dx ' 

9 9 

defined as p z = {x — i)M + ; p M Q , and the free mass Mq acquires the familiar form Mq 



m Q + P 2 + \j m 1 P + V 2 -. with m sp being the mass of the quark-spectator. 
For the sake of brevity we shall present the results for the two representative ET wave 
functions corresponding to the AL1 || and DSR 0] quark models. These functions re- 
sult from solving either the Schrodinger equation (non-relativistic kinematics T = p 2 /2fi re d, 
Hred = mQm sp /(mQ + m sp ); this is the case AL1) or the Salpeter equation (semi-relativistic 
kinematics T = Mq; this is the case DSR). The AL1 potential is the usual Coulomb+linear 
potential supplemented with a gaussian hyperfine term whose range is mass dependent. The 
DSR potential is much more sophisticated one since in addition to the previous terms, it con- 
tains also effects from instantons (which are necessary to describe light pseudoscalar mesons), 
and from finite quark size. The DSR potential is then convoluted with a quark density to take 
into account finite size effects for quarks. The detailed form of this potential is given in Ref. 
[frO| . The parameters of the AL1 and DSR potentials have been determined with help of a fit 
procedure on many experimental meson resonances, from it to B. 

We choose these two different prescriptions for the kinematics to see whether our results 
are sensitive to some "relativization" at the level of the ET wave function. The main difference 
between w(p 2 ) in the AL1 and DSR models relies in the behaviour at high value of the internal 
momentum. The AL1 model results in the soft wave function similar to that of the ISGW2 
model P| with < p 2 >= 0.235 GeV 2 for B° and < p 2 >= 1.075 GeV 2 for B c , while the DSR 



wave function exibits high momentum components leading to < p 2 >= 0.517 GeV 2 for B° and < 

p 2 >= 1.388 GeV 2 for B c . As a result the DSR distribution function is broader than that for the 

case ALL This is illustrated in Fig.l which shows the distribution functions F b (x) for B° and 

i 

B c mesons. Note that F% (x) = F% (1 — x). In terms of the mean value x = J xF(x)dx and the 

o 

i 

variance a 2 = J(x - x) 2 F(x)dx, one obtains x B o = 0.9 (AL1), 0.89 (DSR), cx| = 0.0028 (AL1), 
o 

0.0065 (DSR), x Bc = 0.72 (AL1 and DSR), and a 2 Bc = 0.0062 (AL1), 0.0081 (DSR). The 
constituent quark masses in a relativized quark models are systematically lower that those 
determined using the non-relativistic models. In our case m u = = 0.315, m s = 0.577, m c = 
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1.836, m b = 5.227 GeV (case AL1) and m u = m d = 0.221, m s = 0.434, m c = 1.686, m b = 5.074 
GeV (case DSR). Note that the constituent quark masses m b and m c satisfy approximately the 
relation mj = m c + 3.4 GeV, which is consistent with the well known formula relating the pole 
masses m bjPole and m cpo ; e in Heavy Quark Effective Theory. 

We are now ready to present an overview over different B c decays and their relative impor- 
tance as obtained within the framework we are advocating. To this end we first calculate the 
partial B decay modes corresponding to various underlying quark subprocesses. In Table 1 
we report the B° partial widths for AL1 and DSR LF models. For comparison, we show also 
the results obtained in |/J using the Gaussian-like ansatz of the ISGW2 model which are very 
similar to the AL1 case. 

Our analysis incorporates 54 exclusive SL and NL decays and 29 inclusive decays including 
two baryon-antibaryon channels. The latter ones were calculated using the Stech approach 
[|i~9"| . We have also included the CKM suppressed b — > u contributions with \V u b/V cb \ ~ 0.1. 



The vector and axial form factors for B° — > D(D*)£i>£ and B° —* ir(p)£i>£ transitions have been 



calculated using the formalism of Refs. ||20|| . When calculating the inclusive rates the important 
question is which value to use for the hadron threshold • This value is very badly defined 
theoretically. In our calculations we have adopted two "natural" choices: M% = Mp + m w 
and M^ = My + m^. In Table 1 we show only the results obtained using M$° = My + 
For MP = M P + we obtain 10% increasing of the calculated B° rate. For each case 
we then fix the effective value of \V c b\ by the requirement that the measured B° meson lifetime 
rgo ~ 1.56 ps is obtained. Our goal here is not to establish a new value of \V cb \, but rather 
to illustrate how our approach works. Moreover, imposing the \V cb \ constraint strongly reduces 
the dependence of the predicted value of t Bc from the uncertainty related to the choice of the 
continuous cc threshold 0. 

We apply the same strategy to calculate different partial rates and the lifetime of the B c . 
The critical point with regards to this issue is inclusive charm decay. Here the energy release 
is not as comfortably large as it is in the case of bottom decay. As a result, our estimations 
of the inclusive charm decay should be more sensitive to a hadronization model. However, this 
decay contributes only « 10% to the total c — > s rate. For this reason we did not include any 
hadronization corrections in our calculations. 

The results for the partial B c decay modes corresponding to the various underlying quark 
subprocesses are collected in Table 2 for M% = My + m n . We also include non-spectator 
contributions from weak annihilation (WA) and Pauli interference (PI). The contribution of 
the annihilation channel is 

r.-E>^(^) , fe) , ( 1 -(^) , -, > (3) 

where we take fs c ~ 0.5 GeV, c T = 1 for the r + z/ r channel and c = (2c + (2/i re(i ) + c_(2/i re(i ) 2 /3 
for the cs channel, with jl red = m b m c / {m b + m c ) being the reduced mass of the be system. 
We use c + {jl re d) = 0.8, c_(ji re d) = 1.5. For the T PI we use the expression given in pl| . Note 
that because of the substantial cancellation of weak annihilation rate and the effect of the 
Pauli interference diagrams, no significant uncertainty on the lifetime arises from the limited 
knowledge of the decay constant fs c - We find T a = 0.268 ps -1 , T PI = —0.142 ps -1 . 
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The values of other CKM parameters entering calculations are \V SC \ — 0.974, \Vd c \ — \V US \ = 0.221. 
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We also show in Table 2 the results obtained in |7[] for the ISGW2 LF model f\ and in ||, 
using Operator Product Expansion (OPE). Viewing the latter comparison with due caution, 
regarding the model dependence and other uncertainties in the estimation of the decay modes 
as well as the quark mass uncertainty for the inclusive prediction, it is reassuring that the order 
of magnitude comes out to be consistent. Our bound state corrections are numerically larger 
than very small effects found in ||, especially for the c —>■ s transitions. One possible reason 
being is that an analysis of c — > s decays using OPE is of limited validity, since due to smaller 
energy release the convergence of OPE is slower than for b — > c decays. 

For the sum of b — » c spectator contributions we obtain r^ c ) = 0.501 ps^ 1 (AL1) and 
0.557 ps -1 (DSR), while the total c-decay contribution is r^ 5 ) = 0.900 ps -1 (AL1) and 
1.079 ps -1 (DSR). One observes a dominance of the charm decay modes over 6-quark decays. 
Various branching fractions can be also inferred from Table 2. For instance the semileptonic 
branching ratio BR{B C — > EugX) is found to be 9.3%. The absolute rates are T(B C — > £h>iX c5 ) = 
4.0 (4.4) - 10 13 | V c6 | 2 sec" 1 and T(B C -> iv t X sh ) = 8.0 (9.5) • 10 13 1 | 2 sec" 1 for AL1 (DSR) cases, 
respectively. 

Putting everything together one finds T(B C ) = T l ^(B c ) + T C ^ S (B C ) +r a + r P/ = (0.65 ps)" 1 
(AL1), and (0.57 ps)" 1 (DSR). For = M P + m w we obtain T Bc = (0.61 ps)" 1 (AL1) and 
(0.53 ps) -1 (DSR). We consider the dispersion in predicted values of r Bc as a rough measure of 
our theoretical uncertainty in calculation of inclusive decay rates arising both from the model 
dependence of ip(x,p\) and the different choices of the continuum threshold M% • Averaging 
our predictions we obtain 

t Bc = (0.59 ± 0.06) ps. (9) 

This result compares favourably with estimates obtained using OPE || and also agrees with 
the most recent CDF measurement within one standard deviation. 

Finally, we note that the experimental extraction of B c signal in the hadronic background 
requires the reliable estimation of the branching fraction BR(B C — > J/ip + X), because J/i[) 
can be easily identified by its fi + [i~ decay mode, while the experimental registration of the 
final states containing the r\ c or B^> is impeded by the large hadron background. We obtain 
BR(B C -> J/V> + tt) = 0.1% (AL1 and DSR), and BR(B C -> J/ip + X) = 14.8% (AL1), 
13.6% (DSR). For the exclusive B c — > J/ipfi + i / ^ channel whose signature would be quite clean 
experimentally when J /ip decays into a pair of muons, providing a three muon vertex the 
calculated branching fractions are 1.68% and 1.57% for cases AL1 and DSR, respectively. 

In conclusion, adopting a LF constituent quark model we have investigated the partial 
widths and the lifetime of the B c meson. The hadronic form factors and the distribution 
function were calculated from meson wave functions derived from an effective qq interaction 
intended originally to describe the meson mass spectra. In this way the link between B c physics 
and the "spectroscopic" quark models was explicitly established. For numerical estimates we 
have employed the LF functions related to the ET eigenfunctions of the AL1 and DSR quark 
models. In several important aspects our analysis goes beyond the quark model estimations 
obtained previously. In addition to the SL and the two-meson NL exclusive decay modes we 
have included a number of the B c — > H 5c (Hi s )qiq2 exclusive and B c — > X' 5c (X^ s )qiq2 inclusive 
channels. These channels have considerable impact on the predicted overall b —>■ c rate. To 
sum up, the LF constituent quark model makes clear predictions on the global pattern: (i) a 
short B c lifetime ~ 0.6 ps and (ii) a predominance of charm over beauty decays. 

5 The result of Ref. Q corresponds to fg c = 0.42 GeV. After correcting for the value of Jb c used in this 
paper, one gets the result quoted in Table 2. 
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Table 1 . The partial widths of the B° meson (in the units |14b/0.039| 2 ps _1 ) calculated adopt- 
ing the different choices of the LF wave functions. Each SL width is the sum of 2 exclusive 
and one inclusive widths. For a NL decay the width is the sum of 12 exclusive and 6 inclusive 
channels corresponding to the external and internal W decays. The CKM matrix element \V c b\ 
is calculated using the experimental value of r ex . p (£>°) = 0.641 ps _1 . Also shown are the SL 
branching ratios BRsl and the charm counting n c . 



Model 


AL1 


DSR 


ISGW2 


b —> c + ev e 


0.076 


0.072 


0.074 


b^c + fiv^ 


0.075 


0.071 


0.074 


b —> C + TV T 


0.016 


0.016 


0.015 


b —> c + ud 


0.322 


0.311 


0.313 


b — > c + cs 


0.122 


0.137 


0.123 


b —> c + us 


0.016 


0.016 


0.016 


b —> c + cd 


0.006 


0.007 


0.006 


b — > u 


0.007 


0.007 


0.007 


B° -> NA CJ A C Z C 


0.023 


0.031 


0.023 




0.662 


0.668 


0.652 


Vbc 


0.0383 


0.0382 


0.0386 


n c 


1.193 


1.205 


1.198 


BRsl 


11.42% 


10.78% 


11.35% 
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Table 2 . Inclusive partial rates of B c (in units ps l ) for the choice of the continuum threshold 
MP = M v + m*. The b -> c rates are calculated using the effective value of |V^,| from 
Table 1. The c — > s rates are calculated using |V^| = 0.974. T a and r PJ are calculated using 
f Bo = 0.5 GeV. 

Table 2 



Model 


AL1 


DSR 


ISGW2 M 

Li— L 


OPE H 

LI— — U 


6 — > c + ez/ e 


0.058 


0.060 


0.061 


0.075 


6 — > c + uz/» 

r p 


0.058 


0.060 


0.061 


0.075 


6 — > c + rz/ r 


0.013 


0.014 


0.013 


0.018 


6 — > c + ud 


0.244 


0.261 


0.259 


0.310 


b —> c + cs 


0.093 


0.117 


0.102 


0.137 


b —> c + ixs 


0.012 


0.013 


0.013 




6 — ► c + c<i 


0.005 


0.006 


0.006 






0.005 


0.007 


0.005 






0.013 


0.019 


0.011 






0.501 


0.557 


0.531 


0.615 


c — > 5 + ez/ e 


0.078 


0.090 


0.081 


0.162 


C S + /i^ 


0.073 


0.085 


0.077 


0.162 


C — > 5 + IX J 


0.688 


0.814 


0.698 


0.905 


C — ► S + 1£S 


0.024 


0.029 


0.023 




c —> d 


0.037 


0.061 


0.036 




p C 


0.900 


1.079 


0.915 


1.229 


6c — > TV T 


0.074 


0.074 


0.074 


0.056 


be — > cs 


0.194 


0.194 


0.194 


0.138 


PI 


-0.142 


-0.142 


-0.142 


-0.124 


^tot 


1.527 


1.762 


1.571 


1.914 


T B C 


0.65 


0.57 


0.64 


0.52 
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Figure 1: The distribution functions Fi(x) = J d 2 p±\ipi(x,p 2 L )\ 2 for the B° and B c mesons 
calculated using the AL1 and DSR LF wave functions (this work) and the ISGW2 LF wave 
function (Ref. @). 
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